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Abstract In this work we present the predictions of the "two-infall model" con- 
cerning the evolution of D, ^He and ^He in the solar vicinity, as well as their 
distribution along the Galactic disk. Our results show that, when adopting de- 
tailed yields taking into account the extra-mixing process in low and intermediate 
mass stars, the problem of the overproduction of ^He by the chemical evolution 
models is solved. The predicted distribution of "^He along the disk is also in agree- 
ment with the observations. We also predict the distributions of D/H, D/0 and 
D/N along the disk, in particular D abundances close to the primordial value are 
predicted in the outer regions of the Galaxy. The predicted D/H, D/O and D/N 
abundances in the local interstellar medium are in agreement with the mean val- 
ues observed by the Far Ultraviolet Spectroscopic Explorer mission, although a 
large spread in the D abundance is present in the data. Finally, by means of our 
chemical evolution model, we can constrain the primordial value of the deuterium 
abundance, and we find a value of (D/H)p < 4 ■ lO"'^ which implies flth' > 0.017, 
in agreement with the values from the Cosmic Microwave Background radiation 
analysis. This value in turn implies a primordial *He abundance Yp > 0.244. 
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1. Introduction 

Chemical evolution models are useful tools to derive the primordial abundances of light 
elements, such as D, "^He and ^He, and to give informations on stellar nucleosynthesis. 
It is well known that while the abundances of metals increase in time with galactic evo- 
lution, the abundance of D is always decreasing since this element is only destroyed in 
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stellar nucleosynthesis (for another view sec MuUan & Linsky 1999, who suggest possible 
D ejection from flares into the stellar wind). On the other hand, '^He and ^He should be 
partially produced, and partially biirned into heavier elements. 

In this paper wc show the predictions of the "two-infall model" (Chiappini et al. 2001), 
concerning the evolution of D, "^He and "^He in the solar vicinity and their distribution 
along the galactic disk. 

This model assumes two main infall episodes for the formation of the halo and part of 
the thick disk, and the thin disk, respectively. The timcscalc for the formation of the thin 
disk is much longer than that of the halo, implying that the infalling gas forming the 
thin disk comes not only from the halo but mainly from the intcrgalactic medium. The 
timescale for the formation of the thin disk is assumed to be a function of the galacto- 
centric distance, leading to an inside-out picture for the Galaxy disk formation. 
The two-infall model differs from other models in the literature mainly in two aspects: it 
considers an almost independent evolution between the halo and thin disk components 
(sec also Pagcl & Tautvaisienc 1995), and it assumes a threshold in the star forma- 
tion process (see van der Hulst et al. 1993; Kcnnicutt 1989, 1998; Martin & Kennicutt 
2001). The model well reproduces the majority of observational constraints about the 
abundances of heavy elements both locally and in the whole disk (Chiappini et al. 1997, 
2001). 

The main novelty of this paper concerns the ^He nucleosynthesis prescriptions. Starting 
from the paper of Iben (1967), stellar models have predicted that low- mass stars are 
strong producers of "^He. However, when stellar production of "^He is included in models 
for the chemical evolution of the Galaxy, no agreement between observed and predicted 
■^Hc abundances can be found. In particular, the models predict an overproduction of 
^He (sec Rood et al. 1976; Dearborn et al. 1996; Prantzos 1996; Chiappini et al. 1997). 
On the other hand, there is evidence for ''He production in stars from the observations of 
this isotope in Galactic Planetary nebulae, such as NGC 3242, (Balser et al. 1997; Balser 
et al. 1999), and also from the observations of the ^^C/^^C ratio, which is related to the 
possible occurrence of "extra-mixing" , with consequent "^He destruction, in low and inter- 
mediate mass stars, the main producers of ^He in the universe. The extra-mixing process 
has been suggested to solve the ''He overproduction problem. Recent stellar models sug- 
gested that thermal instabilities might occur during the red giant branch (RGB) phase 
in low mass stars and during the early asymptotic giant branch (AGE) phase in interme- 
diate mass stars, leading to a mixing between ^He-enriched material from the external 
envelope and the material which suffer H-shell burning. These processes should induce 
^He destruction in favor of ''Li production (Sackmann & Boothroyd 1999a; Palacios et 
al. 2001). Extra-mixing should also decrease the ^^C/^^C ratio (see Sect. 3 for a detailed 
description). This justifies the low observed ^^C/^^C ratio in RGB stars and in Li-rich 
giants. 
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However, not every star in the low and intermediate mass range should suffer from extra- 
mixing. For example, recent Hubble Space Telescope observations of NGC 3242, (see 
Palla et al. 2002), suggests a lower limit ^^C/^'^C > 38, in agreement with standard 
stellar models. Therefore, the lack of the ^"'C line and the presence of the "'He line in 
the spectrum of this planetary nebula suggest that the progenitor star did not undergo 
a phase of deep extra-mixing during the latest stages of its evolution. Charbonnel & do 
Nascimcnto (1998) collected the observations of the carbon isotopic ratio in field and 
cluster giant stars, and their conclusion was that 93% of evolved stars undergo the extra- 
mixing on the Red Giant Branch, and are expected to destroy, at least partially, their 
^He. Extra-mixing along the RGB is supported by recent observations of metal-poor field 
stars (Gratton et al. 2000) and of stars along the RGB phase in the thick disk of the 
Galaxy (Tautvaisiene et al. 2001). 

In the present work we take into account the detailed prescriptions about "^He produc- 
tion/destruction by Sackmann & Boothroyd (1999a, b) and their dependence upon stellar 
mass and metallicity. The "^He production is strictly correlated with the destruction of 
D and the production of '*He. Therefore, wc consider all these elements. Wc compare the 
predictions for the time-evolution of the abundances of D, "^He and *He and for their 
distribution along the disk, with the available data. The latest data on '^He/H, which are 
observations of the hyperfine transition of "^He^ at 8.665 GHz (3.46cm), are from Bania 
et al. (2002), and show a flat distribution around a value of about 1.5- 10^"'' (with the ex- 
ception of the observation of a larger '^He abundance in a few nebulae such as NGC 3242). 
The latest published data for D/H, D/0, D/N, are from the Far Ultraviolet Spectroscopic 
Explorer mission (Moos et al. 2002) and refer to the local interstellar medium (LISM). 
These observations (Jenkins et al. 1999; Sonneborn et al. 2000; Friedman et al. 2002) 
show a large dispersion in the deuterium abundances beyond a few tenth of kpc in the 
solar vicinity which is not easy to justify since these effects are not seen in the oxygen 
and nitrogen abundances. 

The paper is organized as follows: in Sect. 2 we describe the chemical evolution model, in 
Sect. 3 wc discuss in detail the adopted nucleosynthesis prescriptions for D, ^He and ^He. 
In Sect. 4 the results are presented and discussed and finally, in Sect. 5 some conclusions 
are drawn. 

2. The model 

The model assumes that the Galaxy forms out of two main accretion episodes almost com- 
pletely disentangled. During the first episode, the primordial gas collapses very quickly 
and forms the spheroidal components, halo and bulge; during the second episode, the 
thin-disk forms, mainly by accretion of a significant fraction of matter of primordial 
chemical composition plus traces of halo gas. The disk is built-up in the framework of 
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the insidc-out scenario of Galaxy formation, which ensures the formation of abundance 
gradients along the disk (Larson 1976; Matteucci & Francois 1989). As in Matteucci & 
Francois (1989), the Galactic disk is approximated by several independent rings, 2 kpc 
wide, without exchange of matter between them. The rate of accretion of matter in each 
shell is: 

r1V..(R +\ ^ ^^^^ ^ ^^^^ ^-it-t^a.)/rD ^ (1) 



dt 

where '£i{R,t) is the surface mass density of the infalling material, which is assumed 
to have primordial chemical composition; tmax is the time of maximum gas accretion 
onto the disk, coincident with the end of the halo/ thick-disk phase and set here equal 
to 1 Gyr; and tjj are the timescales for mass accretion onto the halo/thick-disk and 
thin-disk components, respectively. In particular, = 0.8 Gyr and, according to the 
inside-out scenario, to{R) = 1.03(i?/kpc) — 1.27 Gyr (see Chiappini et al. 2001). A(R) 
and B(R) are derived by the condition of reproducing the present total surface mass 
density distribution in the solar vicinity (see Matteucci 2001). 
The SFR adopted here has the same formulation as in Chiappini et al. (1997): 

where v{t) is the efficiency of the star formation process, S(_R, t) is the total mass surface 
density at a given radius R and given time I](i?Q, t) is the total mass surface density 
at the solar position, GgasiRii) is the gas surface density normalized to the present 
time total surface mass density in the disk Y,D{R,tQai), tQai = 14 Gyr is the age of the 
Galaxy, and R0 = 8 kpc is the solar galactocentric distance (see Reid 1993). The gas 
surface density exponent is k = 1.5, in order to ensure a good fit to the observational 
constraints in the solar vicinity. This value is also in agreement with the observational 
results of Kennicutt (1998), and with N-body simulation results by Gerritsen & Icke 
(1997). 

The efficiency of star formation is set to v — 2 Gyr^^, for the Galactic halo, whereas 
it is I' = 1 Gyr~^ for the disk; this to ensure the best fit to the observational features 
in the solar vicinity. The efficiency v becomes zero when the gas surface density drops 
below a certain critical threshold. We adopt a threshold density {(Jgas)th — 4 Mq pc~^ 
for the Galactic halo, whereas it is {(Jgas)th — 7 Mq pc~^ for the disk (see Chiappini et 
al. 2001). 

The IMF is that of Scalo (1986), assumed to be constant in time and space (for the effects 
of a variable IMF, see Chiappini et al. 2000). 

The yields of metals are from van den Hoek & Groenewegen (1997) for the low and 
intermediate mass stars; from Thielemann et al. (1996) for the nucleosynthesis from 
SNell; from Thielemann et al. (1993) for the nucleosynthesis from SNela; from Jose & 
Hernanz (1998) for the nucleosynthesis from nova outbursts for the elements ^Li, ^^C 
and ^^N. See Sect. 3 for the detailed nucleosynthesis prescriptions of D, ^He and ^He. 
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Table 1. Model parameters common to all the models. 
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0.8 Gyr, timescale for mass accretion 
onto the halo/thick-disk 

1.03(-R/kpc) — 1.27 Gyr, timescale for mass accretion 
onto the thin-disk 

1 Gyr, the time of maximum gas accretion 
onto the thin-disk 

2 Gyr^^, the efficiency of star formation, 
for the halo/thick-disk 

1 Gyr"'^, the efficiency of star formation, 
for the thin-disk 

~ 4 Mq pc~^, the gas surface density threshold 
below which f = for the halo/thick-disk 

~ 7 Mq pc~^, the gas surface density threshold 

below which i/ = for the thin-disk 

1.5 

1.35, if M < 2 Mq 

1.7, if 2 Mq < M < 80 Mq 

3.5 kpc, scale length for the exponential distribution 

of the disk present surface mass density 

54 Mq pc~^, the present surface mass density of the disk 

8 kpc, solar galactocentric distance 

14 Gyr, age of the Galaxy 



Table 2. Baryon to photon ratio r] and flbh^, using the observed value of D/H in high- 

redshift objects as baryometer, and Standard Big Bang Nucleosynthesis. 

For a comparison, ^bh^MB ~ 0-032lQ;QQg, at 95% Confidence Level, see Jaffe et al. 

(2001). 

O'Meara et al. (2001) Levshakov et al. (2002) 

10"^-(D/H) 3.0 ±0.4 3.75 ±0.25 

10-^°-?? 5.6 ±0.5 4.37 - 5.32 

Q.bh^ 0.0205 ± 0.0018 0.016< Q.bh^ <0.020 



The present time surface mass density distribution of the disk is exponential with scale 
length Rd = 3.5 kpc normalized to T,D{RQ,tGai) = 54 Mq pc^^ (see Romano et al. 
2001). The main parameters are presented in Table 1. 
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2.1. The adopted primordial abundances 

In Table 2 wo show some of the most recent estimates of the primordial D abundance. 
These are mainly obtained from observations of Damped Lyman Alpha (DLA) and 
Lyman limit systems, which are supposed to be objects with N(HI)>2-10^° atoms cm~^ 
and N(HI)>3-10^^ atoms em^^ respectively, and where the D abimdance is probably close 
to the primordial value. In Table 2 we also show the cosmological parameters derived 
using D as a baryomcter from Standard Big Bang Nucleosynthesis. Within the Standard 
Big Bang scenario the deuterium abundance is a monotonic, rapidly decreasing function 
of the baryon density rj. Once the Big Bang Nucleosynthesis begins, D is rapidly burned 
into ^He and *He. Moreover, the higher the baryon ratio, the faster the burning and the 
less D survives (Steigman 2000). This strong baryon abundance dependence combined 
with the simplicity of the evolution of D/H in time (being only destroyed inside stars) 
is responsible for the unique role of deuterium as a baryomcter. Given the above, we 
will consider the different observational estimates of the primordial D/H abundance and, 
for each of these measurements, obtain the primordial values of '^He/H and *He/H by 
comparison with the theoretical predictions of Buries et al. (1999). 

In particular, O'Meara et al. (2001) take the weighted mean of D/H values from ab- 
sorbers toward QSOs, PKS 1937-1009, Q1009+2956, HS0105+1619, and find (D/H)=3.0 
±0.4 ■ 10"'\ whereas Pettini & Bowen (2001) take the weighted mean of the D/H values 
from absorbers toward QSOs, HS0105+1619, Q2206-199, Q0347-3819 and find (D/H)= 
2.2 ±0.2 • 10"'\ D^Odorico et al. (2001) find (D/H)= 2.24±0.67 • 10"^ in the absorber 
toward QSO 0347-3819; the latest analysis of this object by Levshakov at al. (2002) gives 
(D/H)=:3.75 ±0.25 • 10^^. We have also run models with a high initial abundance of 
deuterium (D/H « 10""*) as suggested by earlier observations of absorbers (for example 
Songaila et al. 1994; Carsweh et al. 1994; Rugers & Hogan 1996; Bmdes et al. 1999; Tytlcr 
et al. 1999). These observations suggest a high primordial deuterium, but arc uncertain, 
and can be compromised by contaminating H (sec Kirkman et al. 2001). Moreover, mod- 
els that can reproduce the bulk of the observational data in the Galaxy predict only 
a modest D destruction (see also Tosi et al. 1998), thus implying a low primordial D 
abundance. 

The observations of the Proto Solar Cloud (PSC) (Geiss & Gloeckler 1998), and of the 
LISM (Linsky 1998 and Moos et al. 2002), also represent important constraints to the 
deuterium initial abundance, in the sense that the primordial D abundance should be 
higher than the D abundance in the PSC and this in turn should be higher than in the 
LISM. 

In Table 3 we show the initial abundances of D, ^He and ^He adopted in our models. 
We note that the primordial abundance of ^He/H obtained in this way for cases II to IV 
is within the range estimated by Bania et al. (2002). These authors obtained a value of 
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^He/H = 1.1 ± 0.2 X 10"'' for a HII region located at 16.9 kpc from the Galactic center. 
They suggested this value to be close to the primordial one given to the fact that in 
the outer parts of the disk less contamination by stellar activity is expected. Note that 
the primordial values of "^He/H obtained when assmning a high primordial D are not 
in good agreement with the value suggested by Bania et al. (2002). For *He the values 
shown in Table 3 can be compared with the recent values suggested by Izotov & Thuan 
(1998), Viegas et al. (2000) and Gruenwald et al. (2002). It should be noted that the 
primordial Yp estimates still suffer from important systematic errors (see Pagel 2000). 
The primordial abundances of model I are from Chiappini & Matteucci (2000). 



3. Nucleosynthesis prescriptions for the light elements 

3.1. Deuterium 

In the present models, our choice is to assume that deuterium is only destroyed by stellar 
nucleosynthesis (see for example Mazzitelli & Moretti 1980). 

3.2. 3 He 

^He is produced in low and intermediate mass stars during H-burning, but is also de- 
stroyed to form heavier species. Besides the normal nucleosynthesis processes, other mech- 
anisms can destroy ^He. 

In particular, extra-mixing could destroy a significant amount of ^He in low and inter- 
mediate mass stars, although the physics of the driving mechanism is not yet firmly 
established (see below). 



3.2.1. Extra-mixing in low and intermediate mass stars 

■^Hc burning could be affected by extra-mixing processes in low and intermediate mass 
stars. In fact, extra-mixing is probably related both to the mechanism which leads to the 
formation of Li-rich giants, and to the mechanism which could modify the ^^C/^'^C ratio, 
in the same mass range. In standard models, during the first dredge-up, which follows 
the exhaustion of the hydrogen burning in the core, the convcctivc envelope reaches the 
region with freshly synthesized '^He. In the absence of extra-mixing processes, the '^He 
in the envelope survives in the star until it is ejected into the ISM by stellar winds and 
planetary nebulae. 

In particular, in low-mass stars, which ascend the RGB with a degenerate He core, the 
extra-mixing process could help to connect the ^He enriched envelope to the hotter inner 
region near the H-burning shell; ^He could be then burned via the Cameron & Fowler 
(1971) mechanism, to create ^Li via ^He{a,^y Be, which quickly decays into ^Li via 
Be{e^ ^VeY Li. For the freshly synthesized ^Li to survive, the ^Be must be transported 
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rapidly to cooler regions before the decay occurs. Depending on the mixing mechanism 
and its efficiency, which may vary from star to star, the ''Li may or may not reach the 
stellar surface. Models with extra- mixing (see for example Charbonnel 1995; Charbonnel 
et al. 1998, 2000) show a ^^C/^'^C ratio which is below the standard value. 
The evolutionary track of an intermediate-mass star is different because of the faster 
RGB phase, which inhibits the occurrence of extra-mixing on this short period. However, 
extra-mixing could occur, again via a Cameron & Fowler (1971) like process, during the 
early AGB phase, after the exhaustion of helium burning in the core, followed by the 
He-shell burning. 

In stars with initial mass larger than ~ 4 Mq, where the shift of helium burning from 
the core to a shell causes the extinction of the hydrogen shell, deep- mixing could also be 
inhibited, but hot-bottom burning could again lead to ^Li production (see Sackmann 
& Boothroyd 1999a; Charbonnel & Balachandran 2000; Palacios et al. 2001, for further 
discussion and details). 

Sackmann & Bothroyd (1999a) proposed a model, where deep extra-mixing transports 
envelope material into the outer wing of the hydrogen burning shell, where it undergoes 
partial nuclear processing, and then transports the material back out to the envelope. 
Charbonnel & Balachandran (2000) pointed out that in this way Li-rich giants would be 
produced all along the RGB, whereas their data suggest ^Li production in two distinct 
episodes, at the luminosity bump phase of the low mass stars^ and in the early AGB 
phase of low and intermediate mass stars. 

Palacios et al. (2001), showed that under certain conditions a thermal instability can oc- 
cur, triggered by the nuclear burning of the freshly synthesized ""Li in a ^Li burning shell, 
external to the hydrogen burning shell. After the first dredge-up, at the luminosity bump, 
the external wing of the ^Be peak is the first to cross the molecular weight discontinuity. 
When it is connected with the convective envelope, ''Be produced via ^He{a, Be starts 
to diffuse outwards. Depending on the mixing timescale, the Cameron & Fowler (1971) 
mechanism may or may not operate in this region. As evolution proceeds, ^Be diffuses 
outwards and decays in the radiative layer where ''Li is destroyed by proton capture. In 
this thin region the energy released by Li{p, a)a increases significantly. This could lead 
to a thermal instability strong enough to require a convective transport of the energy. 
As a consequence, ^Be could be instantaneously transported in the external convective 
envelope. This results in a sudden increase of the surface ^Li abundance. 



^ In this case the convective envelope reaches the top of the H-burning shell, so that nucle- 
osynthesis occurs at the bottom of the convective envelope. 
^ The RGB luminosity function "bump" phase should be caused by a slower rate of evolution 

for stars in this phase, because of the contact between the H-burning shell and the H-rich, 
previously mixed, zone (Charbonnel et al. 1998, 2000). 
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However, Palacios et al. have not yet published "'He yields and in the present work we 
adopted the '^He nucleosynthesis prescriptions from Sackmann & Boothroyd (1999a,b). 



3.3. ^He 



All the stars eject part of their (old plus new) ^He content in the interstellar medium. The 
prescriptions about ^He nucleosynthesis are taken from van den Hoek & Groenewegen 
(1997) for the low and intermediate mass stars and from Thielemann, et al. (1996) for 
stars more massive than 10 Mq. 

3.4. The production matrix 

Recent chemical evolution models considering the evolution of the abundance of ^He are 
from Dearborn ct al. (1996 - hereafter DST96), Prantzos (1996), GaUi et al. (1997) and 
Palla et al. (2000). The models of Dearborn et al. (1996), tried to solve the problem posed 
by "^He overproduction in many ways. One way is allowing "^He to be reduced to a value 
of about 10^"'' in lower mass stars, with cutoffs of IMq, or 1.6Mq. However, as pointed 
out by the authors, such low-mass stars have long lifetimes, and the reduction in '^He is 
only effective during recent epochs. Another way is to ignore the new ^He production by 
setting the surviving abundance factor g3 = X^f/X^. equal to 1 below 2 Mq, or between 
1 and 2 Mq, where X3. and Xs^ are the initial and the final "^He stellar abundance, 
respectively. However, in this case is difRcult to explain the planetary nebulae data (see 
GaUi et al. 1997). 

Prantzos (1996) assumed the following prescriptions: stars with M > 5Mq are always 
assumed to deplete their initial D + "^He, whereas lower mass stars are assumed to 
maintain their initial D plus '^Hc, plus the "^He produced by the p — p chains. In this 
way, Prantzos (1996) succeeded in avoiding ^He overproduction just by setting to zero 
the p — p contribution. However, these prescriptions do not consider the dependence of 
extra-mixing on stellar mass and metallicity. 

In the present work we take into account extra-mixing, and its dependence from the 
stellar mass, metallicity as well as the "^He content of the proto-stellar cloud. In order 
to account for the ^He nucleosynthesis prescription in our chemical evolution model, we 
adopt the formalism of Talbot & Arnett (1971, 1973). In particular, we define the fraction 
of a star of mass m which is ejected back into the interstellar medium in the form of 
isotopic species i, Run, 



mass. The matrix specifies the mass fraction of the star in which material initially in the 




(3) 



where Qnnj is a time independent "production matrix", and is defined for each stellar 
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form of species j is ejected as species i^, and Xi{t) are the abundance mass fractions. 
Rrrii (t) is a mass fraction, and should be positive, whereas the single elements Qrm j can 
also be negative. It is possible to think about negative Qmi j a-s the consequence of the 
fact that the species i is mainly destroyed to form heavier species. 

According to Talbot & Arnett (1973), we define w2, which is the mass fraction ejected in 
the form of newly created D, and our choice is to set this quantity equal to zero, which 
means that all the D present in the stars is destroyed. The quantity q3 is instead the mass 
fraction within which any original '^He is converted to "^He or heavier species. Talbot & 
Arnett (1973) definition of q3 is q3= max{d,{).QQ) (based on the results of Iben 1967), 
where d is the mass fraction which remains as a stellar remnant. Finally, the quantity 
w3 is the mass fraction ejected in the form of newly created ^He. 

Here we redefine q3 and w3 by including extra-mixing and hot-bottom burning in low 
and intermediate mass stars. 



3.4.1. Redefinition of the elements of the production matrix related to '^He 

^He destruction factors, from Boothroyd (1999)^, have been used to define q3 for stellar 

masses below QMq. 

The new definition of q3 is: 

g3 = 1 - 53, (4) 

where g3 — X^^/X^.. The initial abundance of ^He, X^., is a function of time and it 
depends on the chemical evolution model adopted. In fact, X-^. ~ ^^ig^s + '^1'^ ' -^"^tga, ' 
where -^3^^^^ and -^2^^^, are respectively the interstellar gas abundances of ^He and D 
which are predicted by the model. X^^ is the final stellar abundance of ^He, and it is 
influenced by possible extra-mixing and hot-bottom burning, depending on the stellar 
mass, metallicity and on the initial composition X3 . . 

Attention should be payed to the dependence of g3 from the stellar abundance of ■^He. 
The build-up of the '^He peak during the Main Sequence is independent from the initial 
^He abundance, which is the '^He abundance in the protostellar nebula (following the 
pre main sequence burning). A fraction of the ^He produced during the MS phase is 
then added to the convective envelope after the dredge-up. However, if the initial ''He 
abundance was high, then the added '^He, from the pocket created during the MS phase, 
makes a smaller fractional increase relative to the initial value. Finally, at least in low 
mass stars, extra-mixing on the RGB destroys a fraction of the ''He present in the envelope 
at that time. In this case, the fraction destroyed depends mainly on the star mass and 
metallicity. Therefore: 

HcfiYial ( -^^initial * /"dredge— up ~t- HcmS pocket) * /"extra— mixing (5) 

^ Diagonal elements represent the mass fraction in which species i is unprocessed and ejected. 
* It has been used the CBP-lo case, see http://www.krl . caltech.edu/~aib/ for more details. 
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Table 3. Initial abundances of the models. 
Model 10^-(D/H)p by number lO^-Xop by mass Yp 10^ • (^He/H)p by number lO^-XsHep by mass 



I 2.9 4.4 0.238 0.9 2 

II 2.5 3.8 0.248 1.0 2.3 

III 3.0 4.5 0.247 1.1 2.5 

IV 3.75 5.7 0.246 1.2 2.7 

V 10 15.2 0.240 2.0 4.6 

VI 30 46.3 0.229 3.0 6.9 



where /drcdgo-up is a depletion factor due to dredge-up on the RGB, and /extra-mixing is 
a depletion factor due to extra-mixing (see Sackmann & Boothroyd 1999a, b for further 
details). 

For stellar masses larger than 9M0 we adopted the Talbot & Arnett (1973) q3 definition. 
The quantity w3 has been also redefined, according to Sackmann & Boothroyd (1999a). 
This quantity is a function of g3, which is defined between 0.85 and 9 Mq. For stellar 
masses larger than 9Mq, we define g3 = 1 — g3 (see the appendix for the detailed defi- 
nition). For stars with masses below 2.5 Mq wc adopted w3 = w3cBPBoothroyd ^ given 
percentage / of the stars in this mass range and w3standardDST96 ^^^r the remaining per- 
centage of (1 — /) of the stars. In particular, the quantity / is the percentage of stars in 
the low and intermediate mass range, which suffer from cool bottom processing (GBP), 
and w3standardDST96 is t^c quantity w3 predicted by standard stellar models. For masses 
larger than 2.5 Mq, w3 = w3cBPB„othroyd • 

In Table 4 wc report the various amounts of extra-mixing that we adopted in our mod- 
els. Each prescription for extra-mixing is indicated by a letter (A, B, G, D). Model A is 
essentially the model presented in Ghiappini & Matteucci (2000), except for the fact that 
wc are now adopting different yields for massive stars and a primordial helium abun- 
dance of 0.238 (Gruenwald et al. 2002) instead of 0.241 (Viegas et al. 2000). Model B 
assumes extra-mixing in 70% of the stars with M< 2.5 Mq, as suggested by Galli et al. 
(1997); Model C assumes extra-mixing in 93% of stars with M< 2.5 Mq, as suggested 
by Charbonnel & Do Nascimento (1998) and their study of the carbon isotopic ratio 
12q^13q gg^jj g^j^jj cluster giants stars; Model D assumes extra-mixing in 99% of the 
stars which is quite an extreme case. The case with 100% extra-mixing is excluded by 
the observations of high ^He abundances in planetary nebulae such as NGC 3242 (Balser 
et al. 1999). 
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Table 4. Nucleosynthesis parameters of the models. 
Model q3 w3 extra-mixing 



A Talbot & Arnett (1973) DST96 93% 

B new new 70% 

C new new 93% 

D new new 99% 



Table 5. Observed and predicted quantities at Rg,Q and t =14 Gyr. The predictions of 
the models do not differ since the difference concerns mainly the ^He production. It has 
been assumed a Galactic radius of 18 kpc, thus an area for the Galactic plane of w 10^ 
pc^. 

Models Observed (see Matteucci 2001) 



SNIa (century"^) 


0.3 


o.eh^ 


SNII (century- 1) 


0.9 


O.Sh^ 


*(R3,0,t„o„) (Mo pc-2 Gyr-i) 


2.6 


2 - 10 


(^g(R-g,e,'tnow) (Mq pc"^) 


7 


6.6 ±2.5 


Cg 1 (Tt (Rg,©,tnow) 


0.13 


0.09 - 0.21 


Cnnf{Rg,Q,tnow){MQ pc"^ Gyr"^) 


1 


0.3 - 1.5 


Novae Outbursts (yr~^) 


25 


20 - 30 



4. Results 

4.1. The Solar Vicinity 

A good chemical evolution model should be in agreement with what is called the mini- 
mum set of observational constraints, among which the most important is the G-dwarf 
metallicity distribution (Chiappini et al. 1997; Kotoneva et al. 2002). In Table 5, the 
predicted and observed present-day quantities are shown. In particular the SN rate, the 
star formation rate, the surface gas density, the ratio between the surface gas density 
and the surface total density, the infall rate and the nova outburst rate. The H, D, ^He 
and ^He abundances at 9.5 Gyr (assuming 4.5 Gyr as the age of the Proto Solar Cloud) 
predicted by type I, II, III, IV models, are shown in Table 6. The depletion factors for 
D, and the enhancement factors for ^He and for ^He, Xt^^JXo, are shown in Table 7. 
Xtg^j is the abundance by mass of an element at the present time. 
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Table 6. Abundances at 9.5 Gyr, by mass. The observations are from (a) Anders & 
Grevesse (1989) and (b) Geiss & Gloeckler (1998). 



Element 


A-I 


C-I 


C-II 


C-III 


C-IV 


D-I 


Observations 


H 


0.73 


0.73 


0.72 


0.72 


0.72 


0.73 


0.71 ± 0.02 (a) 


D-IO'^ 


3.4 


3.4 


2.9 


3.5 


4.3 


3.4 


3.0 ± 0.5 (b) 


^He-10^ 


2.2 


3.0 


3.1 


3.3 


3.6 


2.7 


3.2 ± 0.2 (b) 


^He 


0.252 


0.252 


0.262 


0.261 


0.260 


0.252 


0.274 ± 0.016 (a; 



Table 7. Depletion factors, for D [Xi^/ Xt^^^)^ and enhancement factors for ^He and ''He 
{Xt(,^J Xq) at the solar vicinity. 



Element 


A-I 


C-I 


C-II 


C-III 


C-IV 


D-I 


D 


1.5 


1.5 


1.5 


1.5 


1.5 


1.5 


^He 


2.5 


3.0 


2.7 


2.6 


2.4 


1.7 


^He 


1.1 


1.1 


1.1 


1.1 


1.1 


1.1 



4.1.1. Abundances of Deuterium 

In Fig. 1 (upper panel) we show the evolution of D/H with time, as predicted by our 
models I, II, III and IV. The predictions seem to fit the data well for the PSC and the 
LISM, although the data for the LISM show a large spread. The data are from Linsky 
(1998), from the seven lines of sight observed from the FUSE satellite (Moos et al. 2000, 
2002) and from other lines of sight observed by Copernicus, IMAPS, lUE and HST (see 
Table 4 in Moos et al. 2002 and references therein). There is evidence of spatial variability 
for the D abundances in the ISM beyond 0.1 kpc (see the data from FUSE observations 
of Feige 110, Friedman et al. 2002, (D/H)= 2.14 ± 0.82 • 10"^ IMAPS observations of 5 
Orionis, Jenkins ct al. 1999, (D/H)= 0.74i°;J2 . 10-5|, and 7^ Velorum, Sonneborn et al. 
2000, (D/H)== 2.18;° ig • 10-^). 

Strangely enough, the underabundance of D toward 5 Ori A is not accompanied by an 
overabundance of N or O relative to H, as expected because of the fact that destruction 
of D implies creation of N and O (see Steigman 2002). 

The reason of this dispersion in the data is not yet clear. Perhaps some elements are 
systematically removed from the gas phase as they are incorporated into interstellar dust 
(Savage & Sembach 1996), but this does not seem the case for N (Meyer et al. 1997). 
Another possibility is to call for infalling material, as in the case of this model, and for 

^ Similar results are found by Bertoldi et al. (1999), who detect an infrared transition within 
the electronic ground state of the deuterated hydrogen molecule, HD. 
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various efficiencies of mixing processes witliin tlie ISM (see Moos et al. 2002 and Lemoine 
et al. 1999 for a discussion of the possible reasons for the interstellar D/H variations). 
In Fig. 1 (lower panel) , we show the predictions of models V and VI with high primordial 
D abundance. It is clear from the figure that a high primordial D abundance should be 
ruled out since it predicts a too high D abundance for the PSC and the LISM, even 
taking into account possible spatial variability for the LISM D abundance. Therefore, 
cases V and VI are ruled out (see also Tosi et al. 1998). 

Fig. 2 shows the evolution of D abundance as a function of [Fe/H] at different galactocen- 
tric distances, as predicted by model C-I. It is clear from the figure that the model depletes 
much more D near the Bulge, whereas D should be quite close to the primordial value in 
the outer regions of the Galaxy (see Sect. 4.2). Chengalur et al. (1997) detected a low sig- 
nificance feature of DI implying (D/II)= 3.9± 1.0 ■ 10^'"' at large galactocentric distances. 
This value is close to the primordial value and consistent with the low D/H measured 
towards some of the high redshift absorbers. Three recent measurements of D towards 
the galactic center are from Polehampton et al. (2002), Lubowich et al. (2000) and Jacq 
et al. (1999). Polehampton et al. (2002) observed the giant molecular cloud Sagittarius 
B2, located near the Galactic Centre, in the far-infrared. They obtained (D/H) in the 
range (0.2 — 11) -10"®. Lubowich et al. (2000) reported the detection of D in a molecular 
cloud at only 10 pc from the Galactic center. They estimated a (D/H)= 1.7 ± 0.3 • 10^*", 
around 10 times below the local interstellar value. Essentially the same value is given by 
Jacq et al. (1999). Apart from these measurements in the Galactic bulge, there are no 
more measurements in the inner regions of the Galactic disk. 



4.1.2. Abundances of ^He 



Figures 3 and 4 show the time-evolution of ''He. In Fig. 3 we show models A, B. C 
and D computed with initial abundances of type I (see Table 3). Type A models, with 
93% extra-mixing, old q3 and w3 definition do not fit the PSC data very well (Geiss 
& Gloeckler 1998). Type B models, with 70% extra- mixing, new q3 and w3 definition, 
clearly overproduces ■'He with respect to both PSC and ISM observed abundances. Type 
C models, with 93% extra-mixing, new q3 and w3 definition, show an earlier rise and are 
in better agreement with the PSC data, whereas type D models, with 99% extra-mixing, 
new q3 and w3 definition, are quite flat. Fig. 4 shows how the time-evolution depends 
on the initial abundance when the nucleosynthesis prescriptions are fixed. Case V, with 
(D/H)p = 10 • 10-5 which implies (^He/H)p = 2.0 • 10-^ disagrees not only with the PSC 
and LISM D abundances, but also with those of ^He. The same is true for case VI. All 
models, except type B models, seem to be acceptable, thus indicating that a percentage 
of extra- mixing larger than 70 % should be assumed (see also Galli et al. 1997). 
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4.1.3. Abundances of *^He 

The predicted solar ^He abundance by type II, III and IV models are in good agreement 
with the observations. Type I models instead, predict a slightly lower *IIe abundance. 

In Fig. 5 we show the model predictions for Y versus O/H. We obtain, for model 
C-IV, AY/A{0/H) = 52 for 360 < 106(O/H)< 660. This value is in good agreement 
with the value of 50.2 ± 3.9 recently reported by Maciel (2001), based on a planetary 
nebula sample of 81 objects. We notice that, although planetary nebulae are interesting 
objects to study the helium evolution in the Galaxy, these objects can have contaminated 
helium abundances. Maciel (2001) took into account the He contamination by planetary 
nebulae progenitor stars and obtained a lower value than the previous AY/A{0/H) and 
AY/AZ estimates based on the same kind of objects (e.g. Chiappini & Maciel 1994). 

Our C-IV model predicts AY/AZ ~ 1.2 if the entire metallicity range is considered, 
whereas it is ~ 2.75 for 0.016 <Z< 0.018. For 0.008 <Z< 0.018, AY/AZ ~ 1.5. The 
AY/AZ ratio is usually determined from photoionized nebulae such as HII regions and 
HII galaxies (e.g. Lequeux et al. 1979; Pagel et al. 1992; Pagel 2000). More recent results 
are from Esteban et al. (1999) who obtained AY/AZ = 2.87 ± 1.0, 2.81 ± 1.0 and 1.87 ± 
0.5 for Orion, M8 and M17 respectively. These authors adopted Yp = 0.240 ± 0.005. This 
estimate is consistent with that of Pagel & Portinari (1998), who found AY/AZ = 3±2, 
from the analysis of the metallicity dependent location of the main sequence of nearby 
stars from Hipparcos data for stars with Z < Zq. Finally, Maciel (2001) estimates 2.0 
< AY/AZ < 2.7 for Yp = 0.24 and 2.9 < AY/AZ < 3.6 for Yp = 0.23. The predictions 
of model C-IV are in good agreement with the above observations. 

4.2. The Galactic Disk 

In this section wc show our results for the predicted radial abundance distributions of D. 
■'Hc and ^He. In Figures 6 and 7 we show the predicted gradients of D/H, D/0 and D/N 
(the local values of these ratios in the ISM are also shown). The predicted gradient for 
D/H is positive and steep, due to the fact that D is destroyed in stars (see also Prantzos 
1996). The evolution of the inner disk regions are faster than the outer parts, which are 
still in the process of forming thus having an almost primordial composition (Fig. 6). In 
Fig. 6 (upper panel) we show the three measurements of D in the Galactic center which 
show a factor of ~ 10 depletion with respect to the solar vicinity value. We predict a 
depletion factor of ~ 5.7 at 2 kpc from the Galactic center and ~ 2 at a distance of 4 
kpc. For large galactocentric distances (> 14 kpc) a value close to the primordial one 
is expected. In fact, D is probably the most sensitive element to radial variations in the 
timescale of disk formation. The gradients for D/0 and D/N are even steeper, because 
when D is destroyed by stellar nucleosynthesis, O and N are produced as a consequence 
(Fig. 7). The predicted gradients for D/H, D/0 and D/N are c;iog(D/H)/dR ~ 0.02 
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dcx/kpc, dlog(D/0)/dR ~ 0.13 dex/kpc, dlog(D/N)/dR c:i 0.14 dex/kpc, respectively, 
over the 4-18 kpc range. 

In Figures 8 and 9 the predicted ''He abundance distribution along the disk is shown. For 
this element, we predict a negative gradient, since it is more produced than destroyed. In 
the inner regions, which are older in the inside-out scenario, there has been more stellar 
activity, and the contribution of low and intermediate mass stars to the '^He enrichment 
of the ISM has been more important than in the outer regions. Fig. 8 shows how the 
■'He distribution is dependent from the percentage of extra-mixing, whereas Fig. 9 shows 
how the distribution is dependent from the adopted "primordial" ''He abundance, once 
the extra-mixing percentage is fixed. It can be seen that the assumption that 99% of the 
low and intermediate mass stars suffer extra-mixing leads to a flatter distribution, with 
-0.02 < (ilog(^He/H)/(iR< -0.01 dex/kpc. The case with 93% extra-mixing predicts 
—0.04 < (ilog(^He/H)/rfR< —0.03 dex/kpc, in agreement with the data from Galactic 
HII regions (Bania et al. 2002). 

Finally, in Fig. 10 the *He distribution along the disk is shown, with a predicted small 
gradient (i(^He/H) /dR « —0.002 kpc~^ over the 4-18 kpc range. The shaded region shows 
the locus of the planetary nebulae data (see Maciel 2001 for the details). Planetary neb- 
ulae data extend up to ~ 12 kpc and show an essentially flat distribution along the disk. 
Taking all of these results into account, we conclude that our best model is C-IV, which 
assumes an extra-mixing percentage of 93%, which is also the percentage suggested by ob- 
servations (Charbonnel & do Nascimento 1998). This model also assumes the primordial 
D/H abundance of 3.75 • 10"'^ which implies (3He/H)p = 1.2 • 10"'"^ and Ip = 0.246. The 
primordial abundance of deuterium is from Levshakov et al. (2002). These authors mea- 
sured one of the highest value for the deuterium abundance within the set of "low-D/H" 
observations of high redshift absorbers. 

5. Conclusions 

In this paper we have calculated the galactic evolution of the abundances of D, ^He, and 
^He. In particular, we modelled the evolution of the ^He abundance taking into account 
in detail the extra-mixing process in low and intermediate mass stars, and its dependence 
on stellar mass and metallicity. We have predicted the time-evolution of the abundances 
of D, ^He, and '^He, and their distribution along the Galactic disk, and studied how these 
abundances depend upon the "primordial" abundances as well as upon the percentage 
of extra-mixing adopted in the models. We have compared these predictions with the 
most recent observations, including the first published data from the FUSE mission and 
selected the best combination of parameters. 

The adopted chemical evolution model, the "two-infall" one, differs from other models 
present in the literature because it assumes that the galactic thin-disk formed out of 
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an infall episode completely disentangled from the formation of the galactic halo, with 
almost no contribution from the pre-enriched halo gas especially in the inner regions. As a 
consequence of this, our predicted deuterium depletion factors are lower than in previous 
models (e.g. Prantzos 1996), especially in the region close to the bulge. The reason for 
this resides in the fact that after the first halo phase, where some D depletion already 
occurs, the formation of the thin-disk implies mostly primordial gas and therefore the D 
abundance in the gas increases again and then declines (see Tosi ct al. 1998). For what 
concerns the predicted abundances of '^He our model differs from previous ones because 
it includes new nucleosynthesis prescription for this element. 
The main conclusions are the following: 

- The redefinition of the production/destruction of '^He in low and intermediate mass 
stars has important consequences on the evolution of the "'He abundance and on its 
distribution along the Galactic disk. In particular, the new models with 93% of extra- 
mixing, show a good agreement with the PSC observations as well as with the observed 
gradient along the disk. We predict a ^He gradient of —0.04 < (ilog(^He/H)/(iR< 
—0.03 dex/kpc in the galactocentric distance range 4-18 kpc. 

- The predicted time evolution of D in the solar vicinity is well reproduced although 
the LISM values measured by FUSE show a large spread. Since D is only destroyed 
during galactic evolution, we predict a positive gradient of this element along the 
disk, in particular (ilog(D/H)/dR~ 0.02 dex/kpc, dlog(D/0)/(iR~ 0.13 dex/kpc, 
(ilog(D/N)/(iR2± 0.14 dex/kpc, in the galactocentric distance range 4-18 kpc. Our 
predicted D/H gradient is less steep than the one obtained by Prantzos (1996) by 
means of a model with a much lower infall rate than the "two-infall" model. This 
is a consequence of the fact that D abundance is very sensitive to the infall/star 
formation history in the Galaxy. It should be noted that the infall/star formation 
history assumed here is requested to fit the majority of the observational constraints 
(Chiappini et al. 1997, 2001). We predict that observations towards the Galactic 
anticenter should provide a D abundance close to the primordial value. We predict a 
depletion factor for D in the solar vicinity of ~ 1.5, whereas at 2 kpc from the Galactic 
center we expect a factor of ~ 5. However, the model of the Galactic disk is probably 
not appropriate to describe the Galactic bulge which should have suffered a much 
higher star formation than we can deduce from the present models (see Matteucci et 
al. 1999). 

- Chemical evolution models can constrain the primordial value of the deuterium abun- 
dance and a value of (D/H)p < 4 • 10~^ is suggested by the present model. This upper 
limit implies, according to Standard Big Bang Nucleosynthesis (Buries et al. 1999), 
rife/i^ > 0.017, in agreement with estimates from the Cosmic Microwave Background 
radiation analysis, which implies flbh^yi^ = 0.032to.oo8J 95% Confidence Level 
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(Jaffe et al. 2001). This limit seems to rule out high D abimdance deteetion in DLA 
absorbers, and implies Yp > 0.244 from SBBN. This is at variance with the latest 
estimate of primordial ^He from Gruenwald et al. (2002), who suggest that the pri- 
mordial ''He should be Yp = 0.238±0.003 (see also Sauer k Jedamzik 2002). Finally, 
we predict an enrichment of ^He relative to metallicity of AY/AZ = 1.5 in the range 
0.008 < Z ^ 0.018, and a gradient along the galactic disk (ilog(''He/H)/dR « -0.007 
dex/kpc, or d(^He/H)/dR —0.002 kpc~^, in the galactocentric distance range 4-18 
kpc. 
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Appendix A: detailed definition of the quantity w3. 

The quantity w3, the mass fraction ejected in the form of newly created ^He, as defined 

from DST96 is: 

if 0.65 Mq < < 2.5 M©, 

ka = (0.00007 + 0.00135- AC^-^) 
kh = (0.55 - 0.3 • ZogioM*) 



^3,* — ^3,i 

w3 = 



2 

A:a + fcb(^3.*- 2.1 -10-4) 



M* • 0.77 

if 2.5 Mq < KU < 25 Mq, 

ka = (0.000111 +^^:^?^- 0.5- 10-6. 

kb = (0.485 - 0.022 -eT^) 

kc = 1.7 • 10"'^e-^.5.io-fi^ 

fca + fcb(^3,*- 2.1- 10-4) _fc^ 



w3 



Ah ■ 0.77 

if 25 Mq < M,, 

ka ^ (0.00011 - 4 • lO-'^M,) 

kb ^ 0.332 

fca + fcb(X3., -2.1-10-4) 
^ M-0J7 ' 

where M* is the stellar mass, Xa^* is the initial stellar '^He abundance by mass (after the 

pre- main-sequence burning of D into ^He - see for example Mazzitelli & Morctti 1980). 

The quantity w3 as defined by Boothroyd (1999), over the entire mass range, is a function 

of g3, which is defined as g3 — X^^/X^. between 0.85 and 9 Mq, and is a function of 

„, M,, and the stellar metallicity Z. 

If 0.85 Mq < < 4 Mq, 

X3 , • g3{X3 M^, Z) ■ (0.95M, - 0.5) 

^3 = — — '■ ; 

M* 

if 4 Mq < M* < 12 Mq, 

• 53(^3,*, M„ Z) ■ (0.9M, - 0.3) 
wd = — ■ ■ : 

M, 

if 12 Mq < M, < 12 Mq, 

^3,* • 33(X3,,, M,, Z) ■ (0.982M, - 1.22) 



w2i 



M* 

For stellar masses greater than 9Mq, we define g3 as 1 — g3, where q3 is the mass fraction 
within which any originalj ^He is converted to ^He or heavier species. 



That is, within the protosteUaj nebula. 
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Figure 1. Predicted and observed (at Scr) abundances of D as function of time (Gyr): 
Proto Solar Cloud data arc from Geiss & Gloeckler (1998), and LISM data are from 
FUSE mission observations along seven lines of sight (see Moos et al. 2002); observations 
from other lines of sight (see table 4 from Moos et al. 2002 and references therein) are 
also shown. The LISM observations are quite sparsed and there is evidence for spatial 
variability of the D abundances (see Jenkins et al. 1999; Sonneborn et al. 2000; Friedman 
et al. 2002). 

Model predictions arc labeled as follows: upper panel C-I (solid), (D/H)p = 2.9 • 10~^; 
C-II (dotted), (D/H)p = 2.5 • 10"^; C-III (short-dashed), (D/H)p = 3.0 • 10"^ C-IV 
(long-dashed), (D/H)p = 3.75- 10"^ lower panel - C-V (solid), (D/H)p = 10- IQ-^; C-VI 
(dotted), (D/H)p = 30 • 10"^ 
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Figure 2. Abundances of D as a function of [Fc/H], at different galactoccntric distances: 
R = 4 kpc (long-dashed), 6 kpc (short-dashed), 10 kpc (dotted), 16 kpc (solid). We show 
the predictions of model C-IV which assumes 93% extra-mixing, new q3 and w3 definition 
and (D/H)p = 3.75 ■ 10^^. There is more D depletion near the Bulge. At 16 kpc there is 
almost no depletion and the predicted value at 14 Gyr is close to the primordial value. 
The predictions from models A, B and D do not differ from the one shown here since the 
difference concerns mainly the ^He nucleosynthesis. 
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Figure 3. Abundances of "'He as a function of time (Gyr) as predicted by models assum- 
ing different amounts of extra- mixing: model A-I (solid), with 93% extra- mixing, old q3 
and w3 definition; model B-1 (dotted), with 70% extra- mixing, new q3 and w3 definition; 
model C-I (short-dashed), with 93% extra- mixing, new q3 and w3 definition; model D-I 
(long-dashed), with 99% extra- mixing, new q3 and w3 definition. 

Also shown are the Proto Solar Cloud and ISM measured abundances (at 3a) from Geiss 
& Gloeckler (1998) and Linsky (f998). The enhancement factors are ~ 2.5 for the model 
A-I, ~ 3.0 for the model C-I and ~ 1.7 for the model D-I. Model B-I clearly overproduces 
^He. Models B, C and D, with redefined q3 and w3, show an earlier rise compared to 
model A. 
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Figure 4. Abundances of ^He as a function of time (Gyr) as predicted by models A, C 
and D. 

Lower panel: predictions of type-D models, with 99% extra-mixing and new q3 and w3 

definition. The curves are labelled as: D-I (sohd), {^Re/R)p = 0.9 ■ 10"'"'; D-II (dotted), 

(3He/H)p = 1.0 • 10-^ D-III (short-dashed), ('^He/H)p = 1.1 • 10^'^ D-IV (long-dashed), 

(3He/H)p = 1.2 ■ 10-^ D-V (dot-short-dashed), (^He/H)p = 2.0 • 10^^ 

Middle panel: predictions of type-C models, with 93% extra-mixing and new q3 and w3 

definition. The curves are labelled as: C-I (solid), (3He/H)p = 0.9 • 10"^ C-II (dotted), 

(3He/H)p = 1.0- 10-^ C-III (short-dashed), (^He/H)p 1.1 • 10"^ C-IV (long-dashed), 

(3He/H)p = 1.2 ■ 10-^ C-V (dot-short-dashed), (^Hc/H)p = 2.0 • 10^^ 

Upper panel: predictions of type- A models, with 93% extra- mixing and old q3 and w3 

definition. The curves are labelled as: A-I (solid), (3He/H)p = 0.9 • 10"^; A-II (dotted), 

(3He/H)p ^ 1.0 • 10-'^ A-III (short-dashed), (3He/H)p = 1.1 • IQ-^ A-IV (long-dashed), 

(3He/H)p = 1.2 ■ 10-5; A-V (dot-short-dashed), (^Hc/H)p = 2.0 ■ 10"'^ 

These plots show the dependence of the time-evolution of the abundance of ^He/H from 

the initial abundance, once the nucleosynthesis prescriptions are fixed. 
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Figure 5. Model C-IV prediction for Y vs 10*^ 0/H. The predicted value for Ar/A(0/H) 
in the indicated range is ~ 52. Maciel (2001) estimated a Ay/A(0/H) = 50.2 ± 3.9 for 
Yp = 0.24 based on a sample of 81 planetary nebulae. 
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Figure 6. Distribution along the disk of D as a function of the Galactoccntric distance, 
at t = 14 Gyr. Models C, with 93% extra-mixing, new q3 and w3 definition are shown. 
Curves are labelled as follows: model C-I (soUd), (D/H)p = 2.9 • IQ-^; C-II (dotted), 
(D/H)p = 2.5 • 10~^ C-IV (short-dashed), (D/H)p = 3.75 ■ 10"^ 

Upper panel: the data point at R = 8 kpc is the weighted mean from FUSE mission 
observations (Moos et al. 2002 at 3(t), along seven lines of sight. The three data points 
in the Galactic center are from Polehampton et al. (2002), Lubowich et al. (2000) and 
Jacq et al. (1999). 

Lower panel: the seven lines of sight from FUSE mission observations. Moos et al. (2002) 
and data, at 3<t, from other lines of sight observed by Copernicus, IMAPS, lUE and HST 
are shown. The predicted gradient for all models is (ilog(D/H)/(iR« 0.02 dex/kpc over 
the 4-18 kpc range. 
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Figure 7. Distribution along the disk of D/^^O and D/^^N, at t = 14 Gyr. Models 
arc labelled as follows: model C-I (solid), (D/H)p = 2.9 • 10"^ C-II (dotted), (D/H)p = 
2.5 • 10-^ C-IV (short-dashed), (D/H)p = 3.75 ■ IQ-^. 

The large symbol represents the weighted mean from FUSE mission observations (Moos 
et al. 2002 at Sa) which probe the LISM within 0.5 kpc. The predicted gradient for 
aU models is dlog(D/0)/dR« 0.13 dex/kpc and cilog(D/N)/dRw O.Udex/kpc over the 
4-18 kpc range. 



C. Chiappini et al.: Evolution of deuterium, ^He and ''He in the Galaxy 



29 



o 



t>0 

o 




10 

R (kpc) 



15 



20 



Figures. Distribution along the disk of ^Rc/R at t = 14 Gyr. Model A-I (solid), with 
93% extra-mixine', old q3 and w3 definition, (3Hc/H)p = 0.9 ■ 10"'^ model C-I (dotted), 
with 93% extra-mixing, new q3 and w3 definition, (■^He/H)p = 0.9 ■ 10^'"'; model D-I 
(short-dashed), with 99% extra-mixing, new q3 and w3 definition, (■'He/H)p = 0.9- 10^^. 
Observations of the hyperfine transition of ■'Hc+ at 8.665 GHz (3.46cm), from Galactic 
HII regions are shown. Abundance data are from Bania et al. (2002) (filled hexagons). 
The shadow area represents their suggested average value of (1.9±0.6) -10"^. 
The predicted gradient is -0.04 < rflog(3He/H)/di? < -0.02 dex/kpc over the 4-18 kpc 
range. The figure shows predictions and data normalized to the mean observed value for 
3He/H in the ISM (Linsky 1998). 
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Figure 9. Predicted ''He/H gradient in the Galactic disk. Upper panel: models C, with 
93% extra-mixing, new q3 and w3 definition C-I (sohd), {^lle/H)p ^ 0.9 • 10"^ C-II 
(dotted), (^He/H)p = 1.0 ■ 10~^; C-IV (short-dashed), (^He/H)p = 1.2 • 10"^ The pre- 
dicted gradient is —0.04 < dlog^He/dR < —0.03 dex/kpc over the 4-18 kpc range. 
Lower panel: models D, with 99% extra-mixing, new q3 and w3 definition - D-I (solid), 
(3Hc/H)p = 0.9 • 10-'^ D-II (dotted), (^He/H)p = 1.0 ■ 10"^ D-IV (short-dashed), 
(■^He/H)p = 1.2 ■ 10^'"'. With the same extra-mixing, the model using q3 and w3 re- 
defined shows a distribution along the disk shifted toward higher values. The predicted 
gradient is -0.02 < dlog('^He/H)/(iR< -0.01 dex/kpc over the 4-18 kpc range. The 
figures show predictions and data normalized to the mean observed value for '^He/H in 
the ISM (Linsky 1998). Data are from Bania et al. (2002) (see Fig.8). 
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Figure 10. ^Hc as a function of the Galactoccntric distance for models C with 93% extra- 
mixing, new q3 and w3 definition: C-I (sohd), Yp = 0.238; C-II (dotted), Yp = 0.248; C-IV 
(short-dashed), Yp = 0.246. The predicted gradient is (i('*Hc/H)/dR« -0.002 dex/kpc 
over the 4-18 kpc range. The shaded region shows the locus of the planetary nebulae 
data of Maciel (2001). 



